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Abstract: We calculate the potential between two static quarks in QCD using 
modified boundary conditions for the perturbative expansion. Through a change of 
the Feynman ie prescription we effectively add a "sea" of gluons to the asymptotic 
states with energies below a given scale A. We find that the standard result for the 
static potential gets corrections of order A^/Q^ both at small and large momentum 
transfers Q^. After resummation of the infrared sensitive corrections we find that 
the running coupling as{Q'^) freezes in the infrared and that the exchanged gluon 
gets an effective tachyonic mass. We verify that identical results are obtained in the 
Coulomb and Feynman gauges. 
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1. Introduction 

In perturbative QCD (PQCD) calculations of S-matrix amplitudes quarks and gluons 
are assumed to form free asymptotic states at the initial and final times, t —>■ ±oo. 
It is recognized that this is at variance with observations - partons actually bind to 
form colour singlet hadrons which are the true asymptotic states. Consequently, the 
applications of perturbation theory are restricted to so-called infrared safe observ- 
ables in processes characterized by a large momentum scale Q. All predictions are 
subject to power corrections (A/Q)", where A ~ 200 MeV is the fundamental QCD 
scale. 

It has been noted [1] that PQCD predictions can nevertheless be successfully 
extrapolated to low scales Q ~ A, assuming that the Q-dependence of the running 
coupling as{Q'^) "freezes" at a hadronic scale of order A. Confinement appears to 
change momentum distributions only in a mild way, with PQCD distributions of 
partons being reflected in those observed for hadrons. This motivates us to study 
whether PQCD can be modified so that its use can be extended to low without 
having to introduce the freezing effects "by hand" . 

Formally, there is considerable freedom in making a perturbative expansion. The 
standard arguments justifying an expansion, namely 

• The initial and final times are taken to infinity along a ray slightly tilted wrt. 
the real axis, and 

• The asymptotic configurations have a non-vanishing overlap with the true 
ground state of the theory 
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allow many choices of in- and out-states. The existence of an overlap with the true 
ground state is in practice an assumption, even in the case of standard PQCD where 
the asymptotic states are taken to be the empty "perturbative vacuum" . Considering 
the central importance of perturbation theory in applications of field theory it seems 
desirable to explore the properties of expansions with different asymptotic states. 

Here we will study the effect of adding gluons to the perturbative vacuum. It 
is natural to consider this since the true QCD ground state is believed to be a 
condensate of gluons. Conceivably, the background gluons may mimic the properties 
of the true gluon condensate sufficiently to make the perturbative expansion express 
some of the confinement physics already at low orders. In any case, the above formal 
arguments justifying such a modified perturbative expansion are as compelling as 
those of standard PQCD. 

The specific modification of the asymptotic state we consider has been called 
the "Perturbative Gluon Condensate" [2] . Background gluons with energies smaller 
than a given scale A are introduced by modifying the Feynman ie prescription of the 
gluon propagator in the following way: 

+ ^ ['^(^o - |k|) + S{ko + |k|)] e(A - |k|) (1.1) 



k^-\-ie {k-\-iey -\- ie 2|k 

e>(|k|-A) , 1 r 1 



A;2 + i£ 2 



{k^-ief-V^ (A;o + i£)2-k^_ 



e(A- |k| 



where l/(/c^ + iz) denotes the ordinary Feynman ie prescription and l/(/c + i^)^ de- 
notes the modified one. As was shown [2] for scalar fields, a perturbative calculation 
of any Green function G using the modified propagator (1.1) is equivalent to a su- 
perposition of standard calculations using Feynman propagators with gluons added 
to the asymptotic states, schematically 

(0|G|0) - f n E (n(^0"M G I Vikg^.r^) ■ (1-2) 

\|k|<A nk=0 / k k 

Here the rik = term corresponds to the unmodified expansion, the c„j^ are known 
constants and the sum is over on-shell gluons of momentum k and energy |k| < A. 
We will show here that gauge invariance is maintained when both gluon and ghost 
propagators are modified according to Eq. (1.1). 

Physically, the modified asymptotic states imply scattering ofi^ the "background" 
gluons which prevents the creation of gluons with |k| < A. Technically this can be 
seen from the sign change (1.1) of ie in the free gluon propagator which removes 
pinches between positive and negative energy poles in loop integrals. For a fermion 
propagator, such a change of ie would be equivalent to filling all fermion (or an- 
tifermion) levels up to a Fermi momentum A, and consequently preventing fermion 
pair production in accordance with the Pauli exclusion principle. We are motivated 
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to study the analogous modification of tlie gluon propagator as a way of avoiding 
the production of soft gluons in perturbation theory. Since we effectively superpose 
calculations with different numbers of background gluons as indicated in Eq. (1.2), 
we need not specify the wave function of such a "Dirac gluon sea" (c/. [3]). We shall 
refer to the physics based on the modified gluon propagator (1.1), with the standard 
Feynman ie prescription for quark propagators, as "Perturbative Gluon Conden- 
sate Dynamics" , or PGCD. Formally, the PGCD expansion appears as justified as 
ordinary PQCD. 

The introduction of a fixed momentum scale A in the PGCD propagator (1.1) 
seems to break Lorentz invariance. The perturbative expansion of the amplitude 
for a given process will depend on the reference frame, since the scale A is frame 
independent. Formally the series sums to the same (Lorentz covariant) result in any 
frame, but the rate of convergence is frame dependent. The situation is in this sense 
analogous to the well-known freedom of choice in the renormalization scale. Physical 
arguments must be used to choose an optimal frame for each process. This is in 
fact commonly done in hadron phenomenology. The non-relativistic quark model 
describes hadrons in their rest frames, whereas the parton model is formulated in 
the infinite momentum frame. 

We should emphasize that the boost properties of bounds states are in general 
extremely complicated [4]. In QED, positronium wave functions and energy levels are 
nearly always evaluated in the rest frame, and most efficiently using non-covariant 
methods such as NRQED [5] . Not even general features such as the Lorentz contrac- 
tion of QED bound states have (to our knowledge) been exphcitly demonstrated in 
perturbation theory. In QCD we face the extra challenge that the gluon condensate 
ground state is boost invariant: the gluons carry momenta of 0{Aqcd) in any frame. 
This feature can clearly not be described using perturbation theory - the best we can 
do is to approximate the true ground state with background gluons whose momenta 
arc the same in any frame, as in Eq. (1.1). 

In this paper we consider the effects of PGCD on the static quark potential. 
This implies an automatic frame choice since the static potential is defined only in 
the "rest frame" of the static sources. We shall not further discuss the important 
and non-trivial question of Lorentz invariance. The question of frame choice for a 
general process is beyond the scope of this paper. 

According to the Kinoshita-Lee-Nauenberg (KLN) theorem [6] all infrared sin- 
gularites cancel if one sums over incoming and outgoing states that are degenerate 
in energy. Our procedure of adding soft gluons to the in- and ott^-states introduces 
a similar smearing of the physical observables. It has in fact been shown [7, 8, 9] 
that the "KLN-cancellations" can be accounted for using a similar modification of 
the ie prescription as the one we study here. As discussed in [9] the effects of the 
KLN-cancellation can be thought of as a "KLN vaccum" and the non- vanishing in- 
teractions with the vacuum as "perturbative condensates" . Thus the physical picture 
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appears similar to the PGCD. The KLN-cancellations are vahd in any field theory 
irrespective of whether there is confinement or not, and the energy- resolution (corre- 
sponding to A) can be arbitrarily small. In our interpretation the scale of soft gluons 
is a physical feature related to the ground state of QCD. 

The purpose of this paper is two-fold. On the one hand we want to investi- 
gate whether the PGCD boundary conditions give a perturbative expansion which 
captures some of the physics of QCD at long distances, while leaving unchanged 
standard perturbative results at short distances. As a first test case we calculate 
the QCD potential between two static colour sources in a colour singlet state [10]. 
We compare the ultraviolet and infrared properties of the static PGCD potential 
with results obtained using ordinary PQCD. The second purpose of this paper is to 
check explicitly that the perturbative gluon condensate framework is gauge invariant. 
Hence we do the calculation both in a physical and in a covariant gauge, namely the 
Coulomb and Feynman gauges. 

2. Calculation of static potential 

The QCD potential V{Q'^) between two static colour sources can be defined in a 
gauge invariant way from a Wilson loop [10]. At lowest order the potential is just 
given by one-gluon exchange, V{Q^) = —Cpg^/Q"^, where is the strong coupling 
and — —Q^ — — is the squared momentum transfer which is purely space-like in 
the static approximation, i.e. qo — 0. The PGCD ie prescription does not change this 
lowest order result since the coupling of the background gluons to a source with large 
mass M is suppressed by A/M. At higher orders the fixed coupling is replaced by 
the running coupling after renormalization. Including all higher-order corrections in 
the running coupling gives an effective charge av{Q^) defined by 



where Cp = {N^ - 1) /2Nc = 4/3 for QCD. In the following we will calculate ay (Q^) 
to one-loop order using the PGCD ie prescription (1.1). For convenience we define 
the one-loop correction II(Q, Qo, A) so that the leading order result is factored out. 



where Qo is the renormalization point, i.e., II{Qq, Qo, A) = 0. 
2.1 Coulomb gauge 

Coulomb gauge is the most natural gauge for calculating the static potential [11], 
although the Feynman rules are not as simple as in a covariant gauge such as Feynman 
gauge. Here we will use the Feynman rules of Coulomb gauge given by Feinberg [12]. 
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Figure 1: One-loop diagrams contributing to the static potential in Coulomb gauge. The 
thick lines represent the static quarks, and the dashed lines the instantaneous Coulomb 
propagators. The curly and double lines represent the A and £^-field propagators, respec- 
tively. There is also a mixed A and £^-field propagator which appears in (c). 



The diagrams contributing to the static potential at one-loop order in Coulomb gauge 
are shown in Fig. 1. For clarity we do not include the contribution from light quarks, 
which is the same as in standard PQCD. 

Using dimensional regularization the contribution to the unrenormalized one- 
loop correction U from the Coulomb self-energy diagram shown in Fig. 1(a) is 



4— n 



q2(q 



k)2 {k + ief 



(2.3) 



where Ca — Nc = 3, n is the number of dimensions (n < 4), // is the arbitrary 
dimensional regularization scale, the subscripts i,j denote the space components 
(i,j = 1,2,3), and the ie prescription for the transverse gluon propagator is given 
in Eq. (1.1). We have written the integrand in 4 dimensions since we will not be 
interested in constant contributions to 11. This corresponds to a specific choice of 
renormalization scheme. The fco-integral is conveniently done in Minkowski space 
using ordinary residue calculus, and vanishes for |k| < A since the poles at /cq = ±|k| 
are then on the same side of the real axis. The result is symmetric under k <-> q — k 
and can be expressed as 




(2.4) 



where the G-functions reflect the modifled ie prescription. 

To simplify the remaining integrations it is convenient to choose 



X — 



|k| 



y 



|q-k| 

Q 



(2.5) 
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as new integration variables, with measure 

rx+l 



/foo rx+l r2TT 

(T-^k^ dx dy dipQ''-^x''-^y 
Jo Ax-M Jo 



(2.6) 



where wc have again dropped terms proportional to {n — 4) in the angular integral. 
The integration over the azimuthal angle </? gives 27r and the remaining integral 
becomes 

4-n 



3C 
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47r2 \QJ 
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/ dxx"-^ / dy 
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-X 



+ 2x^y^ + 2x^ + 2y^-l 



IGx^y^ 



(2.7) 



where X — A/Q. Before evaluating the integral we shall add the contributions from 
the remaining diagrams to the integrand. 

According to the rules given by Feinberg [12], the contribution to H from the 
sum of the vacuum-polarization diagrams in Fig. 1(b) and (c) is 



b+c = 



X 



ir ,,4-n f d'^k ( (q-k),( 

I (2^0^ (^ 

fco^ + |[k^ + (q-k)2] 



(q- k),( q- k)/ 
k)2 _ 



kjkj 



q^{k + i£Y{q — k + isY 
After integrating over ko and (/? and using (2.5) and (2.6) this becomes 



(2.8) 
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(2.9) 



Adding the Coulomb self-energy and vacuum-polarization contributions of Eqs. (2.7) 
and (2.9) gives 



n = 



+ 



g' 



4-n 



n-4 



47r2 \Q J 
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^y^{y^ — x^) 
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(2.10) 



Note that the apparent pole at a; = y cancels between the two terms in the integrand. 
Doing the integrals we find the result for the unrenormalized one-loop correction to 
the static potential, 

g^ [11, 2// 11 1 4 . _ 

In ' + — + -A^ + C 



n(g,/x,A) = Ca 



+ 



6 Q(2A+1) 6 4-n 3 

In 



47r2 ^ 

(2A-l)(4A3 + 2A2-5A + 3). 2A + 1 



12A 



2A-1 



(2.11) 
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Figure 2: One-loop diagrams contributing to the static potential in Feynman gauge. The 

thick lines represent static quark propagators, the curly lines gluon propagators and the 
dotted lines ghost propagators. For clarity the vertices have been marked with dots. 



where X = A/Q and C is a renormalization-scheme-dependent constant. This is the 
main result of our calculation. Before analysing it in more detail we check that we 
get the same result if we do the calculation in Feynman gauge. This will at the 
same time constitute a non-trivial verification of the gauge invariance of the PGCD 
is prescription. 

2.2 Feynman gauge 

The diagrams which contribute to the static potential in Feynman gauge at one- 
loop order are shown in Fig. 2. In addition to the gluon propagator corrections of 
Fig. 2(a-c) there is also the vertex correction of Fig. 2(d), which has a non-Abelian 
contribution that does not cancel against the quark wave-function renormalization, 
as well as the crossed box diagram of Fig. 2(e), which has a non-Abelian part that 
is not part of the iteration of the one-gluon exchange. In a general covariant gauge 
a diagram with a three-gluon vertex also contributes, but it vanishes in Feynman 
gauge. For more details on the diagrams that contribute in Feynman gauge and how 
the iteration of the one-gluon exchange works we refer to Fischler [13]. 

Note that we have included the diagram with a four-gluon-vertex shown in 
Fig. 2(c). In dimensional regularization this diagram does not contribute to the 
logarithmic UV-divergence, only to a quadratic divergence which normally cancels 
against the other two gluon propagator corrections. However, since we are modify- 
ing the ie prescription these cancellations are no longer guaranteed and therefore we 
include all diagrams. 

We again use dimensional regularization and calculate the integrands of all di- 
agrams in 4 dimensions since we are not interested in constant contributions to the 
final expression. The result after performing the numerator and colour algebra is 

ig^ ^ 4_„ f dJ'k + (/,- + (7)-+ + lOA-g 

{k + ieY{k + q + ieY {k + ie)^{k + q + is)'^ 

2 ^4 

, (2.12) 



^ q 



{k + iey {ko + iey{k + iey {ko + iey{k + iey{k + q + ie)' 



7 



where each term corresponds to a specific diagram in Fig. 2. The l/(/co + factor 
in the vertex correction and box diagrams comes from the static quark propagator 
and is not to be confused with the PGCD prescription (1.1). Doing the integrals 
over ko and (f and making the variable substitutions x — |k|/(5 and y — \\s. + q\/Q 
we are left with 



4— n 
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47r2 \Q 

+ 2y^x^ - 3y^ - x'^ + 1 
%^(?/^ — x"^) 



6a;^ + 2a;2|/2 - Sx^ - y^ + 1 
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xQiy - A) 
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(2.13) 



Even though the integrand is different from the one of Eq. (2.10) obtained in Coulomb 
gauge, the final result after the integrals are done only differs from Eq. (2.11) by a 
renormalization-scheme-dependent constant. There is thus full agreement between 
the two calculations. 



3. Discussion of result 



Our renormalized rcsTilt for the one-loop contribution to the static potential using 
the PGCD is prescription is 



n(g,go,A) = Ca 



+ 



TT 

2A-Q 
12A 



11 2A + go^4A2 

"6"^''^ATg" + 3g^ + ^ 

^ A3 A2 A \ 2A + g 



(3.1) 



which is obtained from Eq. (2.11) by making a subtraction a.t Q = Qo- The constant 
C is thus determined by the condition 11(^0, Qo, A) = 0. 

A basic control of the vahdity of Eq. (3.1) is that it agrees with the standard 
PQCD result in the Q ^ oo hmit. For A/Q ^ we get 



n(g,go,A) 



A/O-^O 



TT 



11, Qo 



3Q2 



A2 



+ 




(3.2) 



Thus the ordinary asymptotic freedom [14, 15] result is retained with power-corrections 
A'^/Q^. Returning to the complete expression (3.1) we also note that fl is well defined 
for all finite Q/A, including Q = 2 A. More precisely, fl is continuous at Q = 2A but 
the derivative dfl/dlnQ has an (integrable) singularity at that point. 

The leading power-correction in Eq. (3.2) scales as A^/Q^. By contrast, in the 
operator product expansion one expects a A^/Q^ scaling behaviour (see [16] for a 
phenomenological calculation and [17] for a related discussion). In this sense our 
result is more similar to the gluon propagator in the manifestly gauge dependent 
< > gluon condensate [18], which has been argued recently to have a possible 
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physical meaning [19, 20]. A h? /Q"^ scaling of the power-corrections to the potential 
in momentum space was also found in an infrared renormalon analysis by Beneke [21]. 
In this context we note that it is not possible to make direct comparisons of results 
obtained for large Q with calculations made in coordinate space since the Fourier 
transform from momentum space to coordinate space involves an integral over all 
momenta Q. 

The sign of the power-correction in Eq. (3.2) decreases the running of the cou- 
pling since the sign of (A^/Qq — A^/Q^) is opposite to that of ln((5o/Q)- An opposite 
behaviour, namely infrared sensitive short-distance corrections which lead to a confin- 
ing potential were found recently [22]. Since this calculation was made in coordinate 
space the results cannot be directly compared as explained above. We also note that 
the infrared renormalon analysis cannot predict the sign of the power-correction, 
only its scaling [21]. To see whether the negative sign of the power-correction found 
in Eq. (3.2) gives a freezing coupling or a confining potential we have to study the 
small Q behaviour of Eq. (3.1) since a possible fixed point for the evolution equation 
is at Q = 0. 

Expanding our result (3.1) in the limit Q/A ^ we find 



n(g,go,A) 



Ca 



Qy(Qo) 

TT 



c(go,A) + 2^ + o(g^ 



(3.3) 



where C{Qo,A) is a constant. We note several interesting aspects of this. First of 
all we see that the only infrared-sensitive term is of the form A'^/Q^; all other terms 
are either constant or vanish in the hmit Q/A — > 0. Especially there is no logarith- 
mic Q-dependence in this limit, in other words there is no logarithmic running of 
the coupling for small Q/A. (This can also easily be seen directly from Eq. (3.1).) 
Another interesting property of (3.3) is that the sign of the quadratic infrared di- 
vergence A^/Q^ is opposite to the one found in Eq. (3.2) and thus corresponds to a 
linear confining potential when Fourier-transformed to coordinate space. 

On the other hand, the A'^/Q^ term signals a possible breakdown of our expres- 
sion for the static potential at small Q^. A closer analysis of its origin in the Feynman 
gauge calculation shows that it arises in the diagrams with insertions in the single 
gluon propagator shown in Fig. 2(a-c). Power counting shows that this is also true in 
a general covariant gauge. Since these insertions can be iterated the corresponding 
corrections should be resummed as a geometric series. 



l + U{Q,Qo,A) + -- 
1 + II{Q,Qo,A) + 2Ca 
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aviQl) A^' 



+ 



9 



= [i + ^) + • • •] (3.4) 

where z/^ = 2Ca(^v{Q'q)^^ 1'^ ^ tachyonic effective gluon mass squared, = —v^ 
and n is the remainder of 11 after subtracting the quadratically divergent contribution 
/Q"^. At higher orders in there will be other contributions oc Is? jQ^ which will 
make the effective mass scale dependent. We note that according to Chetyrkin, 
Narison and Zakharov [17] the phenomenology of a tachyonic gluon mass is quite 
successful and suggests ~ 0.5 GeV^. More generally, the tachyonic pole indicates 
a qualitative change with decreasing in the physics described by PGCD. The 
implications of this are beyond the scope of the present paper and require further 
study. 

The remaining one-loop correction 11 can be absorbed into a modified running 
coupling q;v((5^, A^), allowing our result to be expressed as 



V{Q^) = -47rC^ ^ ■ (3.5) 



Since ft goes to a constant as Q/A — > the modified coupling ay (Q^, A^) freezes in 
the infrared. On the other hand, at large Q/A, n agrees with the standard PQCD 
result for II up to power corrections of 0{K^/Q'^). Thus OiviQ^^^^^) equals the 
ordinary ayiQ"^) for large Q/A. 

To see in more detail how ay{Q^ ^ K^) freezes in the infrared it is useful to consider 
the one-loop /3-function for this coupling, 

CtlnC^ TT 

Taking the derivative of 11 with respect to In Q we find 

A-^ A 1 Q 



/3o(A/g) = Ca 





Q AAJ |2A-g| 



(3.7) 



which is plotted in Fig. 3. The figure shows that the running of the coupling has 
essentially ceased for Q ^ A. From this it follows that if z/^ is small compared to A^ 
then the coupling freezes in the infrared before the pole at = is reached. The 
figure also illustrates the logarithmic singularity of /3o at Q = 2A. 



4. Summary and conclusions 

We have explored the freedom to modify the boundary conditions of the perturba- 
tive expansion in QCD. More precisely we considered a specific modification, called 
Perturbative Gluon Condensate Dynamics or PGCD, where a low-energy "sea" of 
gluons is added to the asymptotic states by modifying the is prescription for gluon 
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Figure 3: The one- loop coefHcient (}q{K/Q) of the /^-function for the modified running 
couphng av{Q'^,h?). 

(and ghost) propagators. As a consequence the gluon degrees of freedom freeze below 
a scale A, analogously to the behaviour of fermions in a Landau liquid. The gluon 
sea will scatter high-energy quarks and gluons, preventing them from forming free 
asymptotic states. 

In order to investigate the physical relevance of the PGCD expansion we calcu- 
lated the one-loop correction it to the QCD potential between a static quark - anti- 
quark pair. For large we found that ft is unchanged up to power-corrections of 
0{K^ /Q^). Thus the short distance structure of PGCD agrees with standard PQCD. 
At small Q^, on the other hand, we found infrared-sensitive contributions to ft of 
0{v^/Q^) which after resummation give the gluon a tachyonic mass w^g^g = —v^. 
The remaining part of 11 is constant in the limit Q/K^Q and gives an effective 
coupling (5v'(Q^, A^) which freezes for Q < A. 

Our result may be summarized by the expression for the static potential 

_,.C.^^-4.C.|(l.^....). (4., 

By comparison we recall that at a finite quark density, described by modifying the 
is prescription of the quark propagator, Debye screening generates a positive gluon 
mass squared. In coordinate space the jQ^ correction term in (4.1) corresponds to 
a linear confining potential. The physical interpretation of our results for < 
requires further study. 

Our renormalized one- loop correction (3.1) to the static potential has a non- 
trivial dependence on KjQ. The fact that we obtained the same result in two quite 
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different gauges strongly suggests that the PGCD prescription preserves QCD gauge 
invariance order by order in ag. It would be desirable to prove this more generally. 
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